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Abstract 



This work deals with the existence of optimal solution and the maximum principle for 
optimal control problem governed by Navier-Stokes equations with state constraint in 3-D. 
OO ' Strong results in 2-D also are given. 
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o 

t^ : 1 Introduction 

In this paper, we shall study the optimal control problem 

T / p \ pT 






P) Minimize-/ I \'^{y{t,x) - y^{t,x))\A dxdt + h{u{t))dt; 



subject to 



^-iyAy + {y ■ V)y + Vp = Dou{t) + fo{t), inQx (0, T), 

1/(0) = I/O in fi, 
V-2/ = infix (0,r), 

y = on «9fi X (0, T) 



(1.1) 



"h- y{t)eK, vtG(o,r), (1.2) 



where if is a closed convex subset in 

H={y;ye {L\Q)f, V ■ y = 0,y ■ n = on dQ}. (1.3) 

Here fi is a bounded and open subset of M^ with smooth boundary dQ, T > is a given 
constant, z/ > is the viscosity constant, /q G L^(0,T; (L^(r2))^) is a source field, y{x,t) 
is the velocity vector, p stands for the pressure, Dq G L{U; (L^(f2))^), and u G L^(0,T; U), 
where f/ is a Hilbert space. 

The function h : U —^ (— C)0, +oo] is convex and lower semicontinuous, y° G -^^(0, T; H), 
and ^ G L{V, H), where V = ((//q (fi)))^ fl H. Two cases of physical interest are covered by 
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the cost functional of this form, 

(a)^ = dil: the objective is to minimize the energy of control and the difference between state 

function and object function; 

(b)"^ = (i2Vx : the objective in this case is to minimize the energy of control and regularize the 

smoothness of the state function. 

Physically, the cost functional in case (a) means a regulation of turbulent kinetic energy, while 

in case (b) it means a regulation of the square of vorticity. (See [6,7,8] for a discussion on this 

control problem.) 

Let us introduce some functional spaces and some operators to represent the Navier-Stokes 
equation (1.1) as infinite dimensional differential equations. 

Denote by the symbol || ■ || the norm of the space V, which is defined by 



y 



N 
|2 



V / \Vyi\^dx, 



and by the symbol | ■ | the norm of M^ and (L^(r2))^. We endow the space H with the norm of 
(L^(r2))^, and denote by (■, ■) the scalar product of H, (•, ■)(v,v') the paring between V and its 
dual V with the norm || ■ \\v'. Let A e L{V, V) and b : V x V x V -i^ Rhe defined by: 



N 

{Ay, z) = y^ Vvi- Vzidx, Wy,zeV 



and 

N 



b{y, z,w) = y2 / yiDiZjWjdx, \fy, z,w eV 



respectively, where A = £-■ D{A) = {H\n))^ n V. We define B:V ^V'hy 

{B{y),w) = b{y,y,w), Wy^w eV 

Let f{t) = Pfo{t) and D e L{U, H) be given by D = PDq, where P : (^^(fi))^ ^ H is the 
projection on H. Then we may rewrite the optimal control problem (P) as: 

(P) Mini/ \<^^y^t)-y%t))\'+ [ h{u{t))dt; 

^ Jo Jo 



subject to 



with 



y'it) + uAy{t) + By{t) = Du{t) + /(t), 
?/(0) = 2/0, 



(1.4) 



y{t)eK VtG[0,T] (1.5) 

Since /,£)'« G L'^{0,T; H),yo e V, equation (1.4) has a unique solution y G W^''^{0,T; H)n 
L2(0, T; D{A)) when A^ = 2 while in the case A^ = 3, for each u e L'^{0,T;U), there exists 
<T{u) <T such that (1.4) has a unique solution y(-;u) G W^'^{0,T*; H)nL\0,T*; D{A)) 
for all T* < T{u). Here T{u) is given by 

^^""^ ^ ^Ci[\\yor + C-)\\f + Du\\l.io,T;H)? ^^'^^ 
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where Cq is a positive constant independent of yo, u and v (see [3], p. 261, Th.5.10). In order 
to formulate the optimal control problem governed by such system in terms of strong state 
y{-;u), we observe from (1.6) that for each L > 0, there exists T{L) > 0, such that for any 
T* < T{L) and any u G L'^{0,T; U) with ||-Dm||l2(o,T;//) < L, equation (1.4) has a unique 
solution y{-;u) e ^^^'^(0, T*;H)n L'^{0, T*;D{A)). Therefore, the optimal control problem is 
well-posed in the sense of strong solutions if we consider the admissible control set as a bounded 
subset of L^(0, T; U). Another way to formulate the control problem is in the framework of 
weak solutions to equation (1.4), that is y G ^^ = L'^{0, T; V) D C^{0, T; H) n W^'\0, T; V) 
, satisfying (see [3],p.265,Th.5.12), for each ^ eV 

|(7/(t), ^)(v",v) + My, ^) + b{y, y, ^) = {Du + /, ^)iv',v), a.e. t G (0, T) ^^ ^^ 



2/(0) = 2/0. 

where Cw{0, T; H) is the space of weak continuous functions y : [0, T] — t- //. It is known 
that there exists at least a weak solution to equation (1.4) for each u G L^(0,T; U) (see 
[3],p.265,Th5.12). We shall denote ^u, = {{y,u) G ^ x L^{0,T;U);{y,u) solution to 
(2A),y{t)eKyte[0,T]}. 

The main results of this work are about the existence of optimal solution and maximum 
principle for problem (P) in 3-D. In [5,7,1 1], some existence results are given for optimal con- 
trol problems governed by Navier-Stokes equations, wherein the admissible state functions are 
considered as the strong solutions to Navier-Stokes equations while in the present work, we give 
the existence result in the framework of weak solutions to equation (1.4). In [5,10,11], some 
Pontryagin's maximum principle type results are derived for optimal control problems governed 
by Navier-Stokes equations. The main differences between the present work and works men- 
tioned above are as follows. In this paper, we shall give the maximum principle for problem 
(P) with state constraint of pointwise type, i.e. (1.5), and it is not studied in [5,10,11], wherein 
the types of state constraint involved include type of integral, type of two point boundary and 
periodic type. We shall not only consider the stat constraint set K as a closed convex subset of 
H, but also derive the maximum principle when K is a closed convex subset of V. Since the 
state constraint in the second case is stronger, the proof is more precise, and the corresponding 
result is weaker, but physically, it can be applied in some important examples in fluid mechanics 
which will be given in section 4. This is also one advantage of the results derived in this paper 
over those in the mentioned works. 

The outline of this paper is as follows. In section 2, we give and prove the existence of 
the optimal pair for problem (P) by considering the weak admissible pair set ^i^,. In section 
3, we shall formulate the optimal control problem in terms of strong state function, which is 
different from that in section 2, and we get the first order necessary conditions for problem (P) 
with state constraint in two different cases mentioned above respectively. In section 4, we give 
some examples of state constraint covered by the two cases. 

The following hypothesis will be in effected throughout this paper: 
(i) K C H is a closed convex subset with nonempty interior; 

(ii) ^ G L{V; H), D G L{U- H), y' G 1^(0, T; H f] D{^*^)), f G L\0, T; H), yo E V; 
(iii) h : U —^ (— C)0, +oo] is a convex lower semicontinuous function. Moreover, there exist 
a > and C G M such that 

h{u)>a\u\1j + C,\/ueU. (1.8) 

When we study problem (P) in the case that K is a closed convex subset of V, we need assump- 
tion (ii') which is assumption (ii) together with the assumption D G L{U;V). 



We recall some properties of b(y, z, w) here (see details in [3,9]): 
b{y, z, w) = —b{y, w, z), and there exists a positive constant C, such that 



ma 



\b{y,z,w)\ < C\\y\\mA4m2+i\\w\ 

where mi, m2, m^ are positive numbers, satisfying: 

mi+m2 + m3> |, ifm, ^f,V2e {1,2,3} 
mi + 1712 + 1713 > "2", if 3z G {1,2, 3}, m^ = y 

We note also the interpolation inequality: 

||y|U<q|y|lMl2/|ir+i 

where a = m — I E (0,1). Here || ■ \\m denotes the norm of the Sobolev space H^^iVt). 

We give some definition which will be used throughout this paper. 
Definition 1. Given a Banach space E and its dual space E' , we denote by BV{Q, T; E') the 
space of all functionals y : [0, T] — )• E' with bounded variation. For each uj E BV{0, T; E'), 
we define the continuous functional fiui on C([0, T]; E) by 

l^L^{z) = I {z{t),du{t))(E,E'), Vz E C{[0,T]-E) (1.9) 

Jo 

Here {■,-){e,e') denotes the dual product between E and E', and the integral takes in the 
Riemann-Steiljes sense. The measure /^i^j will be denoted by dcu, and if we assume the space E' 
is reflexive, then we have the Lebesgue decomposition 

du{t) = ujadt + dujsit) (1.10) 

where uJa E L^{^, T; E'), Uadt is the absolutely continuous part of measure du, and the func- 
tional Us E BV{0, T] E') is the singular part of dj. In other words, there exists a closed subset 
6 E [0, T] with the Lebesgue measure zero such that dus = 0, on [0, T] \ 6 (see [2], p.51-p.57). 
Denote by M(0, T; E') the dual space of C([0, T]; E), i.e. the space of aU bounded E'- 
valued measures on [0,T], and notice that yu^ E M{0,T; E'). We denote /C = {y G C([0, T]; £;); 
y{t) G K, Vt G [0, T]}, and define the normal cone to /C at y by 

^dy) = {/" e Af (0, T- E')- ij{y - x) > 0, Vx G /C} (1.11) 

2 Existence results 

By admissible pair we mean {y, u) E S^^^ which satisfies equation (1.4) in the weak sense, i.e. 

(1.7). An optimal pair is an admissible pair which minimizes (P). 

Theorem 1. The optimal control problem (P) has at least one optimal pair {y, u). In 2-D, y is 

strong solution to equation (1.4). 

Proof: When N=3, we denote 

F{y,u) = \f \^[y{t)-y\tWdt+ f h{u{t))dt 
^ Jo Jo 

di=mi{l- [ \<^{y{t)-y\t))\'+ [ h{u{t))dt;{y,u)E^^}. 
^ Jo Jo 



Then there exist (?/„, m„) G <^^, such that 

di<F{un,yn)<di + -. (2.1) 

n 

By (1.8) and (2.1), it follows that {un} is bounded in L^(0, T; U). Hence, there exists at least a 
subsequence which again denoted by {«„}, such that 

Un ^ u weakly in L^(0, T; U). (2.2) 

Multiplying equation 



y'n{t) + pAyn{t) + %„(t) = Dun{t) + /(f), 

?/n(0) = 2/0 

by y„, integrating on (0, t), we get that 

|2/n(t)|' + l^ / \\yn\?ds<Ci + C2 / |2/n(s)|'ds, 

^0 Jo 

and it follows by Gronwall's inequality that 

|2 , ., / ||„, ||2, 



This yields that 



yn -^ y weak* in L°°(0, T; if), weakly in L'^{0, T; 1/), 
Ay„ ^ Ay weakly in L2(0, T; V). 



Hence, 



'0 

Finally, we obtain by (2.6) and (2.7) that 



(2.3) 



\yn{s)\' + u \\ynrdt<C. (2.4) 

'o 



(2.5) 



By the properties of the trilinear function b, we have that 

1 3 

\{Byn,w)^v',v)\ < C\yn\^\\ynp\\w\\, 

and it follows that 

[ \Byn\Ut <C [ Wynfdt < C. (2.6) 

Jo Jo 



\^\Ut<C. (2.7) 



^^|weaklyinLl(0,T;n (2.8) 

Byn -^ r] weakly in lI (0, T; V) (2.9) 

To show that (y, u) satisfies equation (1.4), it remains to show that 7]{t) = By{t), a.e. in (0, T). 
By (2.4), (2.7) and Aubin's compactness theorem(See[3], p. 26, Th.1.20), we obtain that 

yn-^y strongly in L\0, T; H), (2.10) 



and it follows that 

Io\(Byn-By,ij){v',v)\ < Jo {\b{yn-y,yn,ip)\ + \y,b{yn-y,'ip)\)dt ^ n) 

^0 asn^ +CX), Vz^ G L2(0,T;r), 

where Y = {ip e C^{Q,); divi/j = 0}. Hence, ri{t) = By{t), a.e. in (0,T). Since h is convex 
and lower semicontinuous, we obtain that 

di < F{y,u) < liminf FK,2/„) < d, (2.12) 

n— >+oo 

We also have that for each t G [0, T], 3t„ G (0, T), such that y{tn) G K, and 

y(tn) -^ y(t) weakly in H. (2.13) 

Since iT is a closed convex subset of H, it's weakly closed, and this yields that y{t) G K, Vt G 
[0, T]. Hence, {y, u) is an optimal pair for problem (P). |j 

Remark 1: As we stated in Section 1, when A^ = 3, if we assume that the admissible con- 
trol set is a bounded subset of L^(0, T; U), then we can consider the strong solution in a local 
time interval (0, T*). By the similar method applied in the proof of Theorem 1, we can get 
the existence resuk, and the optimal state function y G VF^'^(0,T*; H) fl L2(0,r*; D{A)). 
Moreover, the same result follows when the state constraint set K is a closed convex subset of 
V. 



3 The maximum principle 

To get the maximum principle, we need to consider the strong solution of the Navier-Stokes 
equations. As we mentioned in Section 2, when A = 3, we need to consider the problem of 
such case with bounded admissible control set ^ad = {u E ^^(0, T; [/); ||-D^^||l2(o,T;H) ^ L], 
and then we can consider the strong solution to Navier-Stokes equation in (0, T*), where < 
T* = T{L + 6) < T{L). Here 5 > is a fixed constant, and T{L) is given by (2.3), i.e. 

^^^^ ^ ^Cl[\\yo\\' + m\f\\hio,T;H) + L^W ^^'^^ 

Denote ^{h) = {u e L'^{0, T; U); J^ h{u)dt < +oo}. When A = 3, we shall assume that 

^{h) C '^ad- (3.2) 

With this assumption, we can consider the strong solution in [0,T*] in 3-D without control 
constraint which is included in the definition of the function h inexplicitly. We give an example 
here to show that this assumption can be easily fuUfiUed. Let h{u) be 

^ ^ [ +(X), if \\u\\u > r 

We see that if Jp h{u{t))dt < +oo, then ||M(t)||;7 < r,a.e. in [0,T], so ||w||l2(o,T;;7) < C'(r). 
Since in 2-D, the strong solution to equation (1.4) exists on arbitrary time interval (0, T), such 



assumption is unnecessary. We still denote the interval [0, T*] where assumption (3.2) holds by 
[0,T]. 

We need also the following assumption: 
(iv) There are (5, u) e C(0, T; H) x L2(0, T; U) solution to equation 

~z\t) + uA~z{t) + {B\y*{t))yz{t) = B{y*{t)) + Du{t) + f{t), 
5(0) = yo 

such that z{t) G intK, for t in a dense subset of [0, T] . 

Here y*(t) is the optimal state function for the optimal control problem (P). Inasmuch as 

B{y*)eL\0,T-H)), \{B\y*)z, z)\ < '^Wzf + CM', 
we know that equation (3.3) has a solution z G W^'^{0, T; H) f] L^{0, T; D{A)). 

Theorem 2. Let {y*{t),u*{t)) be the optimal pair for the optimal control problem (P). Then 
under assumptions (i)r^(iv), there are p{t) G L°°(0, T; H) and u){t) G BV{0, T; H), such that: 

D*p{t)edh{u*{t)) a.e.[0,T], (3.4) 

Pit) = - [ U{s,t){^*^{y*{t)-y%t)))ds- [ U{s,t)dco{s), (3.5) 



Jt Jt 

and 

{du{t),y*{t)-x{t))>0,yxe}C (3.6) 

'0 

Here D*, ^*, B'{y*{t))*, are the adjoint operators of D, ^ and B'{y*{t)) respectively, where 
B'{y) is the operator defined by 

{B'{y)z, w) = b{y, z, w) + b{z, y,w), W z,w eV 

We recognize in (3.5) the mild form of the dual equation 

p'{t) = uApit) + {B\y*{t)Y)p{t) + ^*^{y*{t) - y\t)) + ^^(t), a.e. in (0, T) 

p{T) = Q ^^ ' 

Theorem 3 below is the analogue of Theorem 2 under the weaker assumption : 

(v) K is a closed convex subset of V , and there are (5, u) G C(0, T; i/) x L^(0, T; U) solution 

to equation (3.3), such that z{t) G intyi^, for t in a dense subset of [0, T]. 

Here intyA' is the interior of K with respect to topology of V . 
Theorem 3. Let {y* (t) , u* (t)) be the solution for optimal control problem (P). Then under 
assumptions (ii'),(iii),(v), there are p{t) G L°°{f),T]V'),uj{t) G BV{f),T]V'), such that (3.4) 
and (3.5) hold, and (3.6) holds in the sense of 

T 

{du{t),y*{t)-x{t))(V',v) > 0,Vx G /C (3.8) 

We define first the approximating cost functions 

Fx{y,u) = l [ {\^{y-y'^)\^ + \u-u*\l)dt + j (Mv) + hx{u)) dt. (3.9) 
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2.0 



where h\ and tpx are the regularizations of h and ip respectively, that is 

II l|2 I |2 

h^{u) = inf{ l'^~^^l'^ + h{v)- V e f/}, ^x{y) = inf{^^ + ¥'(^); ^^H} (3.10) 
Here c/? is the characteristic function of K, which is defined by 

f cx) iix e H\K 



^^^^ - ^ if X G i^. 

The function Lpx is convex, continuous. Gateaux differentiable, and d^px = Vipx = idip)x, 
which is single-valued (see details in [3], p.48, Th.2.9). Denote 

^ = {(y, u) e C([0, T]; H) X L2(0, T; f/); {y{t),u{t)) satisfy equation (1.4)}. 

We prove first 

Lemma 1. There exists at least one optimal pair {yx, ux) for the optimal control problem: 

(Px) Mm{Fx{y,uy,{y,u)e^}. 

and [yx, ux) -^ {y*, u*) strongly in C([0, T]; H) n L'^{0, T; V) x L'^{0, T; U). Moreover, {yx} 
is bounded in C([0, T]; F) n L2(0, T; D{A)). 

Proof: The existence of the optimal pair follows by Theorem 1 and the arguments in remark 1 . 
We shall show the convergence of the optimal pair {yx, ux) in 3-D, and it's easy to prove that 
the same results hold in 2-D by applying the similar method. Since 

dx = Fxiyx,ux) < Fxiy\u*) < F{y\u*) = d, (3.11) 

we have that /^ hx{ux)dt < C, VA > 0. Since 

hiJ^iux)) < hxiux) < h{ux), 

where Jj^ = (1 + Xdh)^^, we know that /^ h{Jx{ux)) < C, VA > 0, and by assumption (3.2), 
we obtain that \\D{Jx{ux))\\l'^{q,T;H) < L,yX > 0. Since 

hxiux) = ^WdhxiuxWu + M^aM) > ^\dhxiux)\' + aWJ^xMlll + C, 
we have that {A||9/ia(ma)P} is bounded in L^(0, T), and it follows that, 

/ \\ux - J^xMWldt = A / A| dhxiux)\\ldt ^ as A ^ 0. 
Jo Jo 

This implies that 3Ao > 0, such that VA > Aq 



\\Dux\\mo,T;H)<L + -. (3.12) 



Multiplying equation 



y'xit) + uAyxit) + Byx{t) = Dux{t) + f{t) 
2/a(0) = 2/0 
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(3.13) 



by y\, integrating on (0, t), it follows by Gronwall's inequality that 

\yxit)? + / \\yx{t)\?dt < C, VA > 0. (3.14) 

Multiplying equation (3.13) by Ay\, integrating on (0, t), with the inequality 



\{Byx,Ayx)\ < C\\yx\\^Ayx\i , 
we obtain that 

\\yx{tW + ^J^ \Ayx{t)\'dt < Co (ibof + l£\f + Dux?dt + 1 £ \\yx{sWds\ 

Here Co is the same constant as that in (3.1). It follows by Gronwall's inequality that 

\\yx{t)r<m 

where 



.WH^^CSkV.Wefo, " 



z/-3t(/)3(0)y ' V '303(0) 

0(O) = Co(|bof + - / \f + Dux\Ht) 



V JO 

By (3.12) and the definition of T, we have that 

3C73 (iboP + l{\f\h,,,r;H, + iL+ ly))' " ^'^'^^^ ' 
Hence 



\\yx{t)f + [ \Ayx{t)\'dt < C{6), VA > Aq. (3.15) 

We mention here that C{5) is a constant dependent on 6, and we shall denote all the constants 
by C in the following without emphasis. By the properties of b, and (3.14), (3.15), we get 

\\Byx\\LHO,T;H) < C\\Ayx\\LHO,T;H) < C, VA > Ao. (3.16) 

This yields that 



Jxy\\LHo,T;H) < C, VA > Ao. (3.17) 

Therefore, on a subsequence convergent to 0, again denoted by A, we have 



Since 



yx{t) ^ yi{t) strongly in C([0, T; H]) n L\0, T; V) 
Ayx{t) -> Ayi{t) weakly in L'^{0, T; H) 
{yx{t))' ^ y[{t) weakly in ^^(0, T; H) 
ux{t) -^ ui{t) weakly in ^^(0, T; U) 

\Byx - Byi\ < C (jyx - yi\^\Ayx\ + \\yx - yi\\^\Ayx - Ayi\^ 



(3.18) 



by (3.15) and (3.18), we also have 

By^it) -^ By,{t) strongly in L\0, T; H) (3.19) 

So {yi{t) , ui{t)) is a solution to equation (2.1). Moreover, since 

^xiyx) = ^\dvx{yx)\' + viJnvx)) > ^\d^xiyx)\^ 

we know that {\\d(pxiyx)\^} is bounded in L^(0,T), and since dipxiyx) = jivx — Jxivx))^ 
where Jxiux) is defined by the inclusion Jxiux) — yx + ^d(p{Jx{yx)) 3 0, we have 

/ \yx-Jl{yx)\dt<\T[ \\d^x{yx)?dt^Qas\^Q, 
Jo Jo 

so yx - Jxiyx) -^ a.e. (0, T). Since Jf{yx) G K, Vt G [0, T], we have that yi{t) G K. VtG 
[0, T]. By (3.18) we know that 

lim J^{ux) = ui weakly in L'^{0,T; U). (3.20) 

A-i-O 

Since the convex function u -^ j^ h{u)dt is lower semicontinuous, we obtain that 

liminf / hx{ux) > / h{ui)dt. (3.21) 

Inasmuch as 



^^0 JO JO 



limini Fx{yx,ux) < \imFxiy\u*) < F{y*,u*), 

A-i-O A->-0 

it follows by (3.18) and (3.21) that 

1 I {\'^{yi{t)-y\t))\^ + \u^{t)-u*{t)\l)dt+ [ h{ui{t))dt<\immiFx{yx,ux) 

2 io Jo ^^0 

<l£ my*{t) - y'm') dt + £ h{u\t))dt 
<\f (I'^lyiW - y\t))?) dt + f h{u,{t))dt. 

This yields that ui = u* , yi = y* and ux{t) -)■ u*{t) strongly in L^(0, T; U). 

Lemma 2. Let zx{t) be the solution to the equation: 

z'^{t) + uAzxit) + {B\yx{t)))zx{t) = B{yx{t)) + Du{t) + /(f), 
z\{Q) = yo- 



(3.22) 



Then Zx{t) — )■ z{t) strongly in C{[0,T]; H) fl L'^{0,T]V), where {z{t),u{t)) is defined in 
equation (3.3), and yx{t) is the the optimal solution in lemma 1. 
Proof: Multiplying equation (3.22) by zx{t),^Q get 

2 ^kA(t)|' + u\\zx{t)\\'' + h{zx{t),yx{t), zxit)) = biyx{t),yx{t), zxit)) + (/ + Du, zx) 
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Integrating on (0, t), since 

3 1 3 1 

\b{zx,yx,zx)\ < C\\zxp\zx\^\\yx\\,\b{yx,yx,zx\ < CH^/aH^ I^j/aI^ I^aI, 
we obtain by Young's inequality that 

kA(t)P + ^ t\\zx{s)rds<C^ f\zx{s)\'ds + C,. 
^ Jo Jo 

It follows by Gronwall's inequality that 

\zxit)\^+ I \\zx{tWdt<C. (3.23) 



Multiplying equation (3.22) by Azx{t), integrating from to t, we get that 

\\\zx{t)r -\\\yor + y j\Azx{s)\'ds 

= I Ky\i.'^)^z\i.'^)^Az\i.'^)) + Kz\{.s),yx{,s),Azx{s)) 
Jo 

+b{yx{s),yxis),Azx{s))ds + (/(s) + Du, zxis))ds. 
Since 

13 11 

\b{yx, zx, Azx) + b{zx,yx, Azx)\ < CdlyAlHI-ZAll^l^^^Al^ + ||2;a||||?/a||2 |AyA|2|A2;A|), 

\b{yx,yx,Azx)\ < C\\yx\\^Ayx\^\Azx\, 
we obtain that 

^lkA(t)f + ^ f \Azx{s)\Hs < C, f\\zx{s)fds + C,. 
^ ^ Jo Jo 

It follows by Gronwall's inequality and (3.23) that 



zx{t)f\ [ \Azxit)\^dt < C (3.24) 

Jo 



Since 

\{iB'iyxit)))zx{t),w)\ < C ( \Azxit)\i + \Ayxit)\i ) \w\ 



we obtain by (3.23) and (3.24) that 

/ \iB'iyxm^M'dt<C (3.25) 

By (3.23), (3.24) and (3.25), we get that 

f \z'^{t)\'^dt<C (3.26) 

Jo 
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Hence 



zx{t) -^ z{t) strongly in C([0, T; H]) n L^{0, T; V) 

Azx{t) -^ Az{t) weakly in L'^{0, T; H) (3.27) 

4(t) ^ z'(t) weakly in L^{0, T; H). 



Since 



\{{B'{y,m^,{t)-iB'iy*it)))-z{t),w)\ 

< Ci\yx - y*\^Azx\ + \\zx - z\\-'\A{zx- z)\-2 ^ (3.28) 

+ \zx - z\2\Ayx\\\zx - zp + \\yx - y*p\A{yx - y*)\^)\w\ 

and yx -^ y*, zx -^ z strongly in L'^{0, T; V) fl C([0, T; H]), we have also 

B'{yx{t)))zx{t) -> {B'{y*{t)))z{t) strongly in L2(0, T; i/) (3.29) 

With above inequalities, passing A to in equation (3.22), we find that z{t) satisfies the equation 
(3.3), and by the uniqueness of the solution, z{t) = 5(t).jj 

We shall denote by U{t, s) and Ux{s, t),0 < s < t < T the evolution operators generated 
by uA + {B'{y*{t))y and uA + {B'{yx{t))y respectively, which are given by U{t, s)^ = ip{t) 
and Ux{s., t)x = ipxit) for < s < t < T, where '0(t) and ilJx{i) are the solutions to 

^\t) = vA^{t) + {B\y*{t))yij{t), 

i^{s) = C 

and 

V^IW = uAMt) + (B'iyxitWMt), .3 3^. 

^a(s) = ^ 

respectively. It is well known and easily seen that such evolution operators exist. Denote by 
f/*(t, s), Ul{s, t) the respective adjoint operators of U{t, s) and Ux{s, t), which are generated 
by uA + B'{y*{t)) and vA + B'{yx{t)) respectively. By the similar method applied in lemma 
2, we can obtain that 

\\Ul{t,s)\\L^H,H)<C 



Proof of Theorem 2: 

step l:(first order necessary condition for approximate problem) Since (yx, ux) minimize the 

functional Fx{y, u), we know that 

and this yields 

(^*^(yA - 2/°), w^a) + (Vhxiux) +UX- u\ u)u + {dipx{yx).wx) = 0, (3.33) 

p 

where wx = limp_5.o ^^~^^ , {y^, ux + pu) G ^, and wxif) is the solution to the equation 

w'^{t) + vAwx{t) + S'(2/A(t))«;A(t) = /^n, wx{Q) = 0. (3.34) 
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Let pxit) be the solution to the backward dual equation 

p',{t) = ,^Ap,{t) + {B'{yx{t)r)px{t) + ^*^(yA(t) - y%t)) + dMvxit)) ., ,., 

PxiT) = ^■'-'^^ 

By (3. 33), (3. 34) and (3.35), we get by calculation that 

{p',{t),wx{t)) + {-Apx{t) - {B'{yx{t)y)px{t),wx{t)) + {ux-u*,u)u = 0. 

Hence 

{-D*pxit) + Vhxiux) + UX- u\ u)u = 0, ^ueU 

Finally, we obtain that 

D*px{t) = S/hx{ux) + ux-u*{t), a.e. te[0,T] (3.36) 

step2: (pass d(pxiy\),P\,dhx{ux) to limit) By assumption (iv) and lemma 2, we know that 
3p > 0, Ai > s.t. zx{t) + ph e K, for t in a dense subset of [0,T], V|/i| = 1, VA > Ai. For A 
fixed, zx{t) is continuous in [0, T], so there exists a partition of [0, T], 

= ti < t2 < ■ ■ ■ < tN-1 <tN = T 

such that \zx{ti) — 2;A(^i-i)| < f, zxiU) + ph e K,'il < i < N. Moreover, since 

Yl / {dMyx{t)),yx{t) - zxiU) - ph)dt >Y. Vx{yx{t)) - i^x{zx{U) + ph)dt > 0, 
we get that 

rT ^ fti 

p / \d^x{yx)\dt <J2 {d^x{yx{t)),yx{t) - zx{U))dt 



{dipxiyxit)),yx{t) - zx{t))dt + ^ / {dipx{yx{t)), zx{t) - zx{U))dt, 
i=i Ju-^ 

and it follows that 

^j \dipx{yx)\dt < J {dipx{yxit)),yx{t)-zx{t))dt 

{p',{t) - uApxit) - {B'{yx{t)y)px{t) - 'rf*^{yx{t) - y'{t)),yx{t) - Zx{t))dt 


T 

(pa, -y'x - i^Ayx - {B\yx))yx + 4 + ^^-^a + {B\yx))zx) + (^*^(?/a - ?/°), 1/a - zx)dt 

= [ {px^Du- Dux) - {'^*'^{yx-y^),yx- zx)dt 
Jo 

T 

{{Vhx{ux) + ux{t) - u\ u{t) - ux{t))u - C^*^(yA - y'),yx - zx)) dt 
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< f [hxiux)-hx{u)]dt + C<C. (3.37) 

Jo 

We set ujx{t) = /J dipx{yx{s))ds, t e [0, T], then by (3.37) and the Helly theorem (see [2],p.58, 
Th.3.5), we know that there exists a function to G BV{[0, T]; H), and a sequence convergent to 
0, again denoted by A, such that 

ujx{t) -^ uj{t) weakly in H for every t E [0, T] (3.38) 

and 



rT pT 



J {dipx{yxis)),xis))ds^ J {duj{s),x{s)),yxeCi[t,T];H),yte[0,T]. (3.39) 
Multiplying equation (3.35) by signpA(^) = |(^' we get that 

:TTbAW = n — r + — ^ — \ — + \ — \ + \ — \ • 

dt \px\ \px\ \px\ \px\ 

Smce\b{px{t),yx{t),px{t))\ < C|pA(t)|5||pA(t)||^||?/A(t)||, we get by Young's inequality that 



1 3 



\Px{t)\ \P\it)\^ 2 \px{t)\ 

Integrating on (t, T), we obtain by Young's inequality that 

|PA(t)| + ^y" \\px{s)\\ds<C^+C2J \px{s)\ds. 

It follows by Gronwall's inequality that 

\\PxmL^iO,T;H)<C, (3.40) 

and so by Alaoglu's theorem, 

Px{t)^p{t) w*-L^{0,T;H). (3.41) 

By (3.32), we infer that 

Uxis,t)'^*'^{yxis)-y%s))ds ^ f U{s,t)^*'^{y*{s)-y%s))ds weakly in H, (3.42) 
t Jt 

and by (3.38), we have 

/ Ux{s,t)d(px{yx{s))ds ^ U{s,t)du{s) weakly in B. (3.43) 

Finally, by (3.41), (3.42) and (3.43), we obtain that 

p{t) = -j U{s,t)i^*^iy*it)-y\t)))ds- J U{s,t)du{s). 
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This means that p(t) satisfies equation (3.5). Since 

T 

{difxiyxit)), yx{t) - x{t))dt > ipxivxit)) - ifxixit)) > 0, Vx G /C, 



by (3.39), we can pass A to to get 

r-T 

{duj{t),y''{t)-x{t))>0, 



so (3.6) holds. To complete the proof, it remains to proof (3.4). By (3.36) and the definition of 
dhx, we have that 

' {hx{ux) - hx{v))dt < [ {D*px- {ux-u*),ux-v)udtyv e L\0,T;U). (3.44) 
Jo 

Remembering that hx{u) < h{u), we obtain by (3.21) and (3.44) that 

f {h{u*) - h{v))dt < I {D*p,u*-v)udtyvEL\0,T;U). (3.45) 

Jo Jo 

This implies the pointwise inequality: 

{D*p, u* - ij)u > h{u*) - h{v), a.e. in (0, T). 
This shows that D*p{t) E dh{u*{t)), a.e. in (0, T) jt 

Proof of Theorem 3: Step 1 is the same as the proof of theorem 2. To pass Pxii), dhx{ux{t)) 

and dipxiyxit)) to limit, we need to prove the following lemma first: 

Lemma 3. zx(t) — >■ z(t) strongly in C([0,T]; V), where zxif) and z{t) are the solutions to 

equation (3.22) and equation (3.3) respectively. 

Proof: We have 

(zxit) - my + ^A{z',{t) - z{t)) + (B'iyximzxit) - {B\y\t))yz{t) = B{yx{t)) - B{y*{t)) 

(3.46) 
Multiplying equation (3.46) by A{zx{t) — z{t)), integrating on (0, t), we get that 

l\\zx{t)-~z{t)f + '^j\A{zx{s)~~z{s))\'ds 

= - b{yxis)-y*{s),zx{s),Azx{s)-Az{s))- b{y*{s), zx{s)-z{s). ^^xis)- Az{s))ds 
Jo Jo 

-b{zx{s) - z{s)iyxis), Azx{s) - Az{s)) - b{z{s), yx{s) - y*{s),Azx{s) - Az{s))ds 

+ / {B{yx{s)) - B{y*{s)), Azx{s) - AS{s))ds 
Jo 

<C (\yx - y*\^^\\Azx\\\A{zx - z)\ + \\zx - zf^lAzx - Az\l^ ds 
+C / (jzx - z\^\zx - zW^Azx - A2\\Ayx\ + \\yx - y*\\^Ayx - Ay*\^A{zx - 5)|) ds 
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+ / {B{y,{s)) - B{y*{s)), Az,(s) - A~z{s))ds 
Jo 

< Ci / \\zx - zW'^ds + C2 ( / \Byx - By*\'^ds + \\yx - y*\\c{io,T];V) + INa - z\\c{[o,T;H]) j 

+^ / \Azx{s) - A~z{s))\'ds 
^ Jo 

Denote C2 (^/J \B{yx{s)) - B{y*{s))\'^ds + \\yx - y*\\c{io,T];V) + I^a - z\\c{[o,T;H])J by ^a- By 
the latter inequality and Gronwall's inequality, it follows that 

sup \\zx{t) - z{t)f < 6xe^'^ ^ as A ^ (3.47) 

t6[0,T] 

Indeed, since B{yx{t)) — > B{y*(t)) strongly in L'^{0,T; H), and zx ^ z strongly in 
C([0, T]; H), it suffices to prove that yx -^ y* strongly in C([0, T];V). We have 

{yx{t) - y*{t)y + uAiyxit) - y*{t)) + Byx{t) - By*{t) = Dux{t) - Du*{t), 

yxiO) - y*iO) = ^^"^^^ 

Multiplying equation (3.48) by Ayx{t) — Ay*(t), integrating on (0, t). It follows that 

l\\yx{t) - y*{t)r + v j' \Ayx{s) - Ay*{s)\'ds 

<Ci f \\yx{t)-y*{t))fds + C2 f {\Dux{s) - Du*{s)\^ + \\yx- y*\\c(io,T;H])) ds 
Jo Jo 

+ ^y \Ayx{s)-Ay*{s)\'ds. 

Since ux{t) -> u*{t) strongly in L^(0, T; H), yx -)■ y* strongly in C([0, T]; H), we obtain by 
Gronwall's inequality that 

\\yx{t)-y*{t)f < C (\\Dux{t) - Du\t)\\l,^^,.^.H^ + \\yx - 2/*||c([o,T;H])) dte''^^ ^ Oas A ^ 0. 

This shows that 

yx -^ y* strongly in C([0, T];V). (3.49) 

So (3.47) holds. We complete the proof of Lemma 3. 

Applying the similar method, we can get the following result 

\\U*x{t,s)Uiv,v)<C 

U*{s,t)^^U*{s,t)C mC{[0,T];V)y^eV ^^'''''^ 

Now we come back to pass px{t),dhx{ux{t)) and dipx{yx{t)) to limit. 

By assumption (v) and Lemma 3, we know that 3p > 0, Aq > s.t. zx{t) + ph E K, 
for t in a dense subset of [0,T], VA < Aq, V||/i|| = 1. By the similar arguments to the proof 
of Theorem 2, we can get that {dipx{yx)} is bounded in L^(0,T; V). We denote ujx{t) = 
Jq dipx{yxis))ds, t E [0, T], then we know by the Helly theorem that there exist a functional 
u E BV{[0, T]; V), and a sequence convergent to 0, again denoted by A, such that 

uxit) -^ u{t) weakly in V for every t E [0,T] (3.51) 
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and 

{dipx{yx{s)),x{s))^v',v)ds^ J {du{s),x{s))^v',v), \/x e C{[t,T];V),\/t e [0,T]. 

(3.52) 
Multiply equation (3.35) by j. — T^ffp- in the sense of the dual product between V and V, denote 

qx{t) = A^^px{t), we have 

-j:\\PxW\\v' = -[] — 77771 ^ H — TTVii 

dt \\Px{t)\\v' \\Px{t)\\v' 

{^*^{y*{t)-y\t)),qx{t)) , {dvx{yx).qx{t)) 



+ 



WpxmW' Wpxmv 



Since 



Kqx.yx.Px) ^ hxl^^e\Px{t)\ 1^ \p,\^ 1 

lu, II ^ ^ ||„ /^N|| <<^\\px\\v' iTir:U<£<-, 

wpxWv \\pxmv> wvxmvf 



Kyx,q\,P\) . ^11 III \px 
< C\\px\\v' 



we obtain by Young's inequality that 

\\p>^mv' + ^ r j^j^ds < c, + c, r\\pxis)\\v'ds+^r-i^j^ds. 

Jt \\PX{S)\\V' Jt ^ Jt \\Px[s)\\v' 

It follows by Gronwall's inequality that 

lbA(t)||L-(o,T;y')<C, (3.53) 

and so by Alaoglu's theorem, 

px{t)^p{t) w*-L^{0,T;V') (3.54) 

By (3.18) and (3.50), we have 

f Uxis,t)^*'>f{yx{s)-y\s))ds ^ f U{s,t)'^*^{y{s)-y^{s))ds weakly inV . (3.55) 
By (3.51), we get 

/ Ux{s,t)d(px{y\{s))ds ^ U{s,t)du{s) weakly in v. (3.56) 

Finally, we obtain by (3.54), (3.55) and (3.56) that 

p{t) = -j Uis,t)i^*^iy*it)-y\t)))ds- J Uis,t)duis). 

This means that p(t) satisfies equation (3.5). Moreover, (3.4) and (3.8) also hold by applying 
the same arguments as that in the proof of theoerm 1. [j 
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We shall consider the reflexive Banach space E as H or V, and denote by (-, ■) the dual 
product between E and it's dual of E (When E = H, it is the scalar product in H), by || ■ || 
the norm of E. Under the hypothesis of Theorem 2 or the hypothesis of Theorem 3, We give a 
corollary here: 

Corollary 1. Let the pair [y*, u*) be the optimal pair in problem (P), then there exist uj{t)E 
BV{[0, T]; E') andp satisfying along with y*, u*, equations (3.4),(3.5),(3.6) (or(3.8)) and 

Uait) G NKiy*it)), a.e.t E (0, T) (3.57) 

du, e ^K.{y*) (3.58) 

HeretJa(t) is the weak derivative of cij(t), and dug is the singularpart of measure rfoj. NK{y*{t)) 
is the normal cone to K at y*{t), and ^K{y*) is the normal cone to /C at y* which is precised in 
definition 1. 
Proof: Let to be arbitrary but fixed in (0,T). For y E K and e > 0, define the function y,. 

y*{t) for|t-to|>e 

y,{t) = <((!- e-\h - t))y + e-\t^ - t)y*{to - e) for t G [to - e, to] 
{1 - e-\t - to))y + s-\t - to)y*{to + e) for tG [to,to + £]. 

Obviously y^ is continuous from [0,T] to E and ^^(t) E K,\ft E [0, T]. By (4.6)(or (4.8)), we 
have 

/ {u{t),y*{t)-y,{t))dt+ [ {dus,y*-ye)>0. (3.59) 

^0 Jo 

We set pe(t) = e~^{y*{t) — yeit)). If to happens to be a Lebesgue point for the function Ua, 
then by an elementary calculation involving the definition of ye, we get 



lim / {uait),pe{t))dt = K(to), ?/*(to) - y). (3.60) 

Inasmuch as y*{t) — ^^(t) = outside [to — e,to + e], we have 

rto+£ 

{dus, y* -ye) = / {dujs, y* - ye)- 

Jto-e 

On the other hand, for each ry > 0, there exists {j/j*^}ili C E, and ctj^ G C([0, T]) such that 

TV 

\\y*{t)-Y,y:,(^rnm<VioitE[Q,T] 

i=l 

We set z^{t) = y*{t) - Xlii y*,ai^(t), then we have 

rto+e 

I / idUs,z^)\<{Vsito + e)-Vsito-e))-snp{\\zJ;\t-to\<e} 

Jto-e 

where Vs{t) is the variation of Us on the interval [0,t]. Since Vs is a.e. differentiable on (0,T), 
we may assume that 



"to+e 
'to-e 
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limsupe:^^ / {dUg, Zjj) < Ct], (3.61) 

e-t-O Jtn-e 



where C is independent of 7]. Now, we have 

rto+e N N pto+e 

JtQ~e i^i i^i Jto 



TV 



i=l 

where V^ir?(0 i^ the variation of (w^, y*^) on interval [0,t] and 7,^ = sup |aj^(t)|. Since the weak 
derivative of Ug is zero a.e. on (0,T). We may infer that: 



-Vi^{t) = 0, a.e. in(0,r), 



and therefore we may assume that 



r-io+e ^ 



Ume-1 /" " (rfa;,, Vy;«,^(t)) = 0, Vr^ > 0. (3.62) 

By (3.61) and (3.62) we have 

lime"^ / {dUs,Pe) = 0. 
^^0 Jo 

By (3.59) and (3.60) we get 

{ua{to),y*{to)-y)>0, a.e.toG(0,T). (3.63) 

Since y is arbitrary, (3.57) holds. 

To conclude the proof it remains to be shown that dug G 'yVj^{y*), that is Jg {dcos, y* — x) > 
0, Vx G J^. Let (9 be the support of the singular measure dus- Then for any e > 0, there 
exists an open subset U of (0,T), s.t. O dU and m(W) < e, where m is the lebesgue measure. 
Let p G Co~(M) be such that < p < l,p = 1 on C» and p = on (0,T) \ W. We set 
y^ = px + (1 — p)?/*, where a; G =J^ is arbitrary. By (3.6)(or(3.8)), we have 

f K(t), y*{t) - y%t))dt + [ {dcOs, y* - y') > 0. (3.64) 

JO Jo 

Since y* - y^ = on (0, T) \ U, we obtain that 

' {ua{t),y*{t)-f)dt\< f \cua\dt<5{e), 
Ju 

where lime_^o S{e) = 0. On the other hand, since dcOs = on (0, T) \0 and p = 1 on C, we 
see that 

{dus, y* - y') = / {dus, y* -x). 



It follows that Jq {dus, y* — x) > —S{e). Since e is arbitrary, (3.58) holds, ft 
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4 Examples 

In this section, we shall give some applications of the above results in some special cases of 
state constraints wherein Theorem 1 and Theorem 2 can be applied. 

Example 1. Let K be the set K = {y E H; J^ \y(x)\'^dx < p^}, then K is a closed convex set 
in H, since 

WmWcmnH) < C{\\B{y*{t)) + Du{t) + fmmo,T;H))) 

so it is feasible to apply theorem 2 to get the necessary condition of the optimal control pair 
after checking whether condition (iv) is satisfied or not. The set K physically gives a constraint 
on the turbulence kinetic energy, which is usually bounded instead of infinite. In this case, the 
maximum principle can be described as following: 

D*p{t)edh{u*{t)) a.e.[0,T] (4.1) 

p[{t) = lyAp^it) + {B\y*{t)y)p,{t)+^*^{y*{t) - y%t))+Ua{t), a.e. in (0,r) 

Pi{T) = 

(4.2) 
p',{t) = uAp^it) + {B'{y*{t)y)p2{t) + dus, a.e. in (0, T) 
_ P2{T) = ^^-'^ 

Moreover, 

^a(t) G NK{y*{t)) = {X{t)y*{t)- X{t) > 0, a.e. in (0, T)} (4.4) 

Here pi,P2 is the decomposition of p, that is p{t) = pi{t) + P2it). Since uja{t) G -^^(0, T; H), 
dujs is the singular part of the measure du, we know that equation (4.2) has a strong solution 
Pi E C{[0, T]; H), while equation (4.3) has only a mild solution p(t) = — /^ U{s, t)dws. 
Example 2. Let K be the so called Enstrophy set 



K = {yEV; /" |V x y\'^dx < p^} 
Jn 



where V x y = curl y{x), and it is true that the norm |V x y\ is equivalent to the norm \\y\\ 
in the space V. In fluid mechanics, the enstrophy S{y) = /^ |V x y\'^dx can be interpreted 
as another type of potential density. More precisely, the quantity directly related to the kinetic 
energy in the flow model that corresponds to dissipation effects in the fluid. It is particularly 
useful in the study of turbulent flows, and is often identified in the study of trusters, as well as 
the flame field. Enstrophy set gives a constraint on the vorcity of the fluid motion. Since 

\\m\\ci[o,Ty,v)<Ci\\Biy*it)) + Du{t) + fit)\\L^o^T-,H))) 

it is feasible to apply theorem 3 to get the necessary condition of the optimal control pair after 
checking whether condition (v) is satisfied or not. In this case, the maximum principle can be 
described by (4.1), (4.2) and (4.3). Moreover, 

oJait) E NKiy*it)) = {Xit)Ay*{ty,X{t) > 0,a.e. in (0,r)} 

Example 3. Let K be the so called Helicity set, 

K = {y eV; / {y, curl yfdx + X \Vy\'^dx < p^] 
Jn Jn 
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where A, p are positive constants. In fluid mechanics, helicity is the extent to which corkscrew- 
like motion occurs. If a parcel of fluid is moving, undergoing solid body motion rotating about 
an axis parallel to the direction of motion, it will have helicity. If the rotation is clockwise 
when viewed from ahead of the body, the helicity will be positive, if counterclockwise, it will 
be negative. Helicity is a useful concept in theoretical descriptions of turbulence. Formally, 
helicity is defined as 

H = {y, curl y)dx 
Jq 

The helicity set plays an important role in fluid mechanics, and in particular, it is an invariant set 

of Euler's equation for incompressible fluids(See [4]). This set gives a constraint on the helicity 

and the smoothness of the velocity field. By the same argument as in Example 2, we know that 

it is feasible to apply theorem 3 to get the necessary condition of the optimal pair when the 

state constrained set is Helicity set, and in this case, the maximum principle can be described 

by (4.1), (4.2), (4.3). Moreover, 

C0a{t) e NK{y*{t)) = {X{t){Ay*{t) + cnT\y*);X{t) > 0,a.e. in (0,T)} 
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